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Abstract 

We propose two ways for determining the Green's matrix for problems ad- 
mitting Hamiltonians that have infinite symmetric tridiagonal (i.e. Jacobi) 
matrix form on some basis representation. In addition to the recurrence re- 
lation comming from the Jacobi-matrix, the first approach also requires the 
matrix elements of the Green's operator between the first elements of the 
basis. In the second approach the recurrence relation is solved directly by 
continued fractions and the solution is continued analytically to the whole 
complex plane. Both approaches are illustrated with the non-trivial but cal- 
culable example of the D-dimensional Coulomb Green's matrix. We give the 
corresponding formulas for the D-dimensional harmonic oscillator as well. 
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I. INTRODUCTION 



Green's operators play a central role in theoretical physics, especially in quantum me- 
chanics, since the fundamental equations are formulated as integral equations containing 
Green's operators in their kernels. Integral equation formalisms have an advantage over 
those based on differential equations because they automatically incorporate the boundary 
conditions. In spite of this fact differential equations are more extensively used in practi- 
cal calculations. The reason certainly is that the Green's operators occurring in integral 
equations are much more complicated than the corresponding terms in the Hamiltonian. 

A possible way of compensating this drawback is using a representation in which the 
Green's operator appears in a simple form. In this respect the momentum-space representa- 
tion is rather appealing, as the free Green's operator is very simple there. This is the main 
reason why momentum-space techniques are so frequently used and also why they are ca- 
pable of coping with complicated integral equations like the Faddeev equations (see Ref. |TJ 
for a review). 

The free Green's operator can also be given analytically between harmonic oscillator 
(HO) states 0. This allowed the construction of a flexible method for solving the Lippmann- 
Schwinger equation in HO-space, which contains the free Green's operator in its kernel ||. 
The representation of operators on a finite subset of a countable basis, such as the HO 
basis, turns the Lippmann-Schwinger equation into a matrix equation. The completeness of 
the basis ensures the convergence of the method. Actually, this is equivalent to a separable 
expansion of the potential (see e.g. Ref. [|J for a review). In this approach only the potential 
term is approximated, the terms in the Green's operator (the kinetic energy term in the case 
of the free Green's operator) are not. Thus, although one is working with finite matrices, 
the solution possesses correct asymptotic behavior. 

To give account of the fact that Coulombic asymptotics are genuinely different from the 
free one, the kernel of the integral equations describing Coulombic systems should include 
Coulombic, rather than free Green's operators. For the two-body Coulomb Green's operator 
there exists a Hilbert-space basis in which its representation is very simple, namely the 
Coulomb-Sturmian (CS) basis. In CS-space the Coulomb Green's operator can be given by 
simple and well-computable special functions, which can be continued analytically to the 
complex plane M. This is also a countable basis, so we have a matrix representation. 

In the past few years a quantum mechanical approximation method for treating Coulomb- 
like interactions in two-body calculations was developed along this line. The analytic 
Coulomb Green's matrix allows the extension of the method to resonant- and scattering-state 
calculations 0. Since only the asymptotically irrelevant short-range interaction is approx- 
imated, the correct (two-body) Coulomb asymptotics is guaranteed. The corresponding 
computer codes for solving two-body bound-, resonant- and scattering-state problems were 
also published |7|]. 

Recently the CS-space representation approach to the Faddeev integral equations has 
been applied to solving the three-body bound- and scattering-state problem with Coulomb 
interactions [§,0. In this formulation of the equations the most crucial point is calculat- 
ing the resolvent of the sum of two independent, thus commuting, two-body Coulombic 
Hamiltonians. This is given by the convolution integral 



2 



(z-h- h)- 1 = ^-i dw(z-w- hi)' 1 (w - h 2 )-\ (1.1) 

Z7T1 JC 

Here the contour C should encircle, in counterclockwise direction, the spectrum of h 2 without 
penetrating into the spectrum of hi . The analytic nature of the two-body Green's matrix 
made the evaluation of the contour integral, also in practice, possible. In fact, the convolution 



integral follows directly from the Dunford-Taylor integral representation of operators [IT 
A function of an operator h is defined as 

f(h) = ^-l dw f(w)(w-h)-\ (1.2) 
2m Jc 

where / should be analytic on C. This way we can calculate complicated functions of 
operators as well. 

Making use of ( |1 . 1|) , we can solve problems which otherwise would amount to solving 
non-separable partial differential equations with unknown boundary conditions. So we be- 
lieve that the analytic representation of simple Green's operators is of extreme importance, 
probably in other fields of physics too, but certainly for the solution of the underlying integral 
equations of quantum mechanics. 

Hamiltonians having Jacobi-matrix forms were also extensively studied in the context of 
the L 2 approach to quantum scattering theory . Here the Hamiltonian is represented 



on an appropriate L 2 basis, which is chosen in such a way, that the asymptotic part of the 
Hamiltonian possesses a Jacobi-matrix form. The resulting three-term recurrence relation 
can be solved analytically, yielding the expansion coefficients of both a "sine-like" S(r) 
and the "cosine-like" solution C(r). The Jacobi-matrix solutions S(r) and C(r) are then 
used to obtain the exact solution to a model scattering problem defined by approximating 
the potential V by its projection V N onto the finite subspace spanned by the first N basis 
functions. So, on the level of physical assumptions the L 2 approach and the methods used in 
Refs. || and j5|-fij] are equivalent in the same sense as the Schrodinger equation is equivalent 
to the Lippmann-Schwinger equation. However, we believe that the approaches starting 
from integral equations, especially those in Refs. are superior since they allow us to 

cope with problems still lacking satisfactory solution, such as three-body scattering problems 
with Coulomb interactions. We note that in the L 2 approach the Green's function of a 
Jacobi-matrix Hamiltonian can also be constructed from the coefficients of the solutions 



S(r) and C(r) [T3|]. However, this construction is applicable only in very exceptional cases. 

In this paper we wish to demonstrate that, if the Hamiltonian appears in a symmetric 
infinite tridiagonal, i.e. Jacobi-matrix form in some basis representation with analytically 
known matrix elements, then the corresponding Green's matrix can be given in terms of 
three-term recurrence relation. We present two independent methods for determining the 
Green's matrix from a three-term recurrence relation. 

In our first method (method A) we consider this relation only as a useful computational 
tool. In addition to the recurrence relation this approach also requires, the matrix element 
of the Green's operator between the first elements of the basis. This often necessitates the 
analytical evaluation of complicated integrals, which restricts its use to exceptional cases 
only. Once this matrix element has been calculated, we can resort to the recurrence relation 
in order to determine the Green's matrix. However, from the numerical point of view the 
recurrence relations can lead to calculational problems and instabilities [IT3 . 
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In our second approach (method B) we propose direct solution of this recurrence relation 
by continued fractions. The richness of the theory of recurrence relations and continued 
fractions enable us to avoid the difficult and strenuous procedure of method A. The inverse 
of the Green's matrix can be evaluated solely from the ratio of two successive elements of 
the recurrence relation. In method B this ratio is provided by a continued fraction. This 
means that, all the above mentioned numerical problems can be avoided, since the recurrence 
relation is completely ignored as a computational tool. 

In Sec. II below we sketch methods A and B. In Sec. Ill we illustrate method A with a 
non-trivial but calculable example, the D-dimensional Coulomb problem. The corresponding 
3-dimensional formulas had been presented earlier 0-0] and were extensively used in both 
two- and three-body calculations. A summary of the relevant mathematical formulas for 
continued fractions and three-term recurrence relations is given in Sec. IV. In Sec. V method 
B is presented using the example of the D-dimensional Coulomb problem. This is followed 
by numerical illustrations in Sec. VI. Finally, in the Appendix the D-dimensional harmonic 
oscillator is considered. We show, that the harmonic oscillator Hamiltonian takes a Jacobi- 
matrix form on a harmonic oscillator basis that has a different frequency parameter. 



II. JACOBI-MATRIX REPRESENTATION 



We first define our representation space for the quantum mechanical problem. Let us 
consider the basis states {\i)} and {\i)}, with i = 0,1,2,..., which form a complete bi- 
orthogonal set, i.e. 

= = Sa (2.1) 



DC 



i = £KX*l = £l*m (2.2) 

i=0 i=0 

Let us start with the defining equation of the Green's operator G corresponding to Hamil- 
tonian H 

1 = (E-H)G, (2.3) 

and apply it to the ket \j) 

\j) = {E-H)G\j). (2.4) 

Inserting a completeness relation between E — H and G, and multiplying form the left by 
the bra (j\ we get 

oo 

% = £(*!(£ -W>(^|j>. (2.5) 

i'=0 

If (i\(E — H)\i') takes a Jacobi-matrix form the infinite sum is reduced only to three terms 
and we arrive at a recurrence relation for the matrix elements Gij = (i\G\j): 



4 



Sij — Ju-iGi-ij + JuGij + Ju + iGi + ij, i — 1, 2, ... ,j — 0, 1, ... , (2-6) 

where J^- = (i\(E — H)\j) are the elements of the Jacobi-matrix. For the % = j = case Eq. 
( 273 ) takes the form 

1 = JoiGiq + Joo^oo- (2-7) 

In method A, if Goo is known, we can calculate Gio from ( |2.7D , and then can continue 
to Gjo using Eq. (|2.6|) . Interchanging the indices in Eq. ( [2.6|) we again get a three-term 
recurrence relation which can be utilized to generate the Gy elements from the known 
Gqj terms. The analytic calculation of Goo together with the application of the recurrence 
relation fl2.6|) with (|2.7|) constitutes the basic idea of method A. 

In method B we compute the Green's matrix without the explicit use of the recurrence 
relation. From the theory of special matrices we know that the inverse of a Jacobi-matrix, 
i.e. the Green's matrix, possesses the property |TH 



G„ I 1 '''''- '][ '. \ I . (2.8) 
Therefore, for j < N and i < N we can write the resolvent equation Q2.5p as 

N 

8ij = /AJii'Gj'i — SiN JiN+l GjN+l) 
i'=0 
N 

= y^AJii' ~ $iN ^i'N JiN+l PN+l/PN)Gi'j 
i'=0 
N 

= y^AJg' ~~ $iN 8i'N JiN+i GoN+i/GoN)Giij, (2.9) 

i'=0 

i.e. the inverse of the truncated N x N G^ Green's matrix is given as 

(Gif ) 1 = Jij — SjN SiN Jnn+i Gqn+i/Gqn- (2.10) 



Eq. (|2.10|) asserts that the inverse of the truncated Green's matrix is determined by the 
elements of the Jacobi-matrix and the ratio Gqn+i/Gqn- This later will be calculated by 
continued fractions derived only from the elements of Jacobi-matrix. This is the basic idea 
of method B. We notice, that in many practical applications, like the solution of Lippmann- 
Schwinger equations, it is directly the (G^P)™ 1 that is really needed (see e.g. Ref. @,[|). 



III. D-DIMENSIONAL COULOMB GREEN'S MATRIX IN 
COULOMB-STURMIAN REPRESENTATION 

Here we define the Coulomb-Sturmian basis and show that on this particular basis the 
D-dimensional Coulomb Hamiltonian has a Jacobi-matrix structure. We also derive the 
corresponding three-term recurrence relation for the Green's matrix and present its solution 
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via method A. In doing so we apply techniques established already for the three-dimensional 
(D = 3) case [|. 

Let us consider the radial Schrodinger equation for the D-dimensional hydrogen atom in 
the /'th partial wave 



H c ip{r) 



Ti 2 ( d 2 1 / D - 3 
2m \ dr 2 r 2 \ 2 



1 + 



D - 1 



Z 

r 



ip(r) = Eip(r). (3.1) 



(See, e.g. Ref. |16| and references.) The bound-state energy spectrum is given by 

mZ 2 

E nr i - 



2h\n r + l + ^) 2 ' 
and the corresponding wave functions are 

r Q T(n r + l) \ 1/2 



(3.2) 



2T(n r + 21 + D — I) 



exp(-%)(a r)^V jL (^- 2 )( aor ) ; 



(3.3) 



where we used the notation a = {{n r + I + -^p-)r ) 1 and r = h 2 /2mZ. 

The Coulomb-Sturmian equation has a structure similar to the eigenvalue equation ( pTT 



1 + 



D - 3 



1 + 



D-l\ (2n + 21 + D - l)6c 



+ 6|U(6 5 ,r) = (3.4) 



and is solved by the Coulomb-Sturmian (CS) functions 



(r\nl) = (j) n i(b s ,r) 



Un + l) 



1/2 



T(n + 2l + D -1) 



exp( 



-bsr)(2bsr) l+I ^ 1 L^ +D - 2 \2bsr), (3.5) 



which are the generalizations of the corresponding CS functions for the three-dimensional 
case P7| . Here bs is a real parameter, thus we have the bound-state CS functions. 

Introducing the notation (r\nl) = <p n i{bs, r ) / r we can express the orthogonality and 
completeness of the CS functions as 



and 



{nl\n'l) = 5 n 



1 = J2 \nl)(nl\ = \nl)(nl\, 



(3.6) 



(3.7) 



ra=0 



n=0 



respectively, confirming that they form a bi-orthonormal basis. 

The overlap of two CS functions can be written in terms of a three-term expression 



(nl\n'l) = (2&5)" 1 [<W(2n + 21 + D - 1) - <W-i((n + l)(n + 2/ + D — 1)) 1/2 
-^n'+i(^(^ + 2/ + J D-2)) 1 / 2 



(3.8) 
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A similar expression holds for the matrix element of H c : 



(nl\H c \n'l) 



4m 



8 nn , [2n + 21 + D - 1 



r bs 



+<W-i((n + l)(n + 21 + D - 1)) 1/2 + <W+i(n(n + 21 + D- 2)) 1/2 j . (3.9) 



^From the above two equations follows the Jacobi-matrix structure of J% n , = (nl\E—H c \n'l). 

The Jacobi-matrix structure immediately implies a three-term recurrence relation (p.6|) 
which, in the method A, can only be solved if Gq is at our disposal. This matrix element 
can be gained by evaluating an integral of the D-dimensional Coulomb-Green's function 

Gf D (r,r',E) = I ^ A<_ i7tf+g _ 1 (-2iAr < )W_ i7>t+# _ 1 (-2ifcr > ), (3.10) 

with the n = n' = CS functions ( p.5|) , where 7 = Z/(2rok) and is the wave number. 
Using the formula by Buchholz 



k) W^ia^M^i^t) 



I f°° ( 1 \ 1 f w\ 2k 

= t(aia 2 ) 2 J exp I — -(ai + a 2 )tcoshwJ J M (t(aia2) 2 sinhw) ( coth — J du> (3.11) 

the integration can be performed analytically |T9|J5[1 and the final result is 



r c 



- iA;) 2 Z + p - l)/2 + i 7 



, , £> - 3 D + 1 

x 2-F1 -/ 7T- + J 7, 1; I + — + 17 + 2; 



'bs + ik 
be - ik 



(3.12) 



IV. CONTINUED FRACTIONS AND THREE-TERM RECURRENCE 

RELATIONS 

Based on the mathematical literature 0,0] we give a brief review of the underlying 
mathematical theorems. Let {a n (z)}f and {& n (z)}^°, a n (z) 7^ 0, be two sequences of com- 
plex valued functions defined on the region D of the complex plane. We define the linear 
fractional transformation s n (w, z) as 

s n {w n , z) = , n > 1, s {w , z) = b (z) + w , (4.1) 

b n (z) + w n 

and 

S n (w n , z) = 5«_i(s n (tt; n , z), z), S (w , z) = s (w , z). (4.2) 
A continued fraction is an ordered pair 
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(({a n (z)},{b n (z)}),{f n (z)}), (4.3) 

where {f n ( z )} is given by 

f n {z)=S n {0,z), n = 0,1,2,3,.... (4.4) 

Here S n (w n , z) is called the nth approximant of the continued fraction with respect to the 
{ w n}^Lo com pl ex series. S n (w n ,z) can be written, using one of the standard notations, as 

Q f v , / v , Oi(z) 2 (z) On(«) K v 

bi{z) + b 2 (z) H h n (z) + U>n 

The convergence of a continued fraction means the convergence of the sequence of approxi- 
mants S n (w n , z) to an extended complex number 

f(z) = \imS n (w n ,z) = b (z) + 7C=i f^|y) , (4-6) 

where 

K oo ( an{z) \ = ai(z) a 2 (z) a n (z) , . 

n=1 \b n (z)J Wiz) + b 2 (z) + ■ ■ ■ + b n (z) + ■ ■ ■■ 1 ' ) 

It should be noted that if f(z) exists for two different sequences of {u n } then f(z) is unique. 
A special class of continued fractions for which the limits 

lim a n (z) = a(z) and lim b n (z) = biz) (4.8) 

n— >oo n— »oo 

exist for all z G D is called limit 1-periodic continued fractions. The fixed points w±(z) of 
the linear fractional transformation 

s(w,z) = lim s n (w n ,z) = a ^ — , (4.9) 

n^oo 5^) + w 

where w = lim^oo w n , are given as the solution of the quadratic equation 

b{z) + w 



w± 



z) = -b(z)/2 ± y/(b(z)/2)* + a(z). (4.11) 



The w±(z) with smaller modulus is called attractive fixed point, while the other one is 
called as repulsive fixed point. Since w±(z) represent the tail of a limit 1-periodic continued 
fraction we can speed up the convergence using the attractive fixed point in the approximant 

S n (w,z). 

The idea of the analytic continuation of the continued fraction f(z) in Eq. (|4.6|) is based 
on the proper choice of {w n } in the approximant S n (w n ,z) [pO~|,pl| . If a continued fraction 



exists in a certain complex region z G D then in many cases it is possible to extend the 
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region of convergence to a larger domain D* D D, where D* depends on the choice of the 
functions w n (z). In the case of limit 1-periodic continued fractions the analytic continuation 
is defined with the help of the fixed points w±(z) of Eq. ( |4.10D as 



f D *(z) = lim S n (w±(z),z). (4.12) 

n— >oo 

The computation of the approximants S n (w±(z), z) might be unstable for certain z, 
which leads to unsatisfactory convergence. This problem can be overcome by using 



the Bauer-Muir transformation [20]. The Bauer-Muir transform of a continued fraction 
b (z) + K (a n (z)/b n (z)) with respect to a sequence of complex numbers {w n }^ =Q is the con- 
tinued fraction d (z) + K (c n (z) / d n (z)) , whose "classical" approximants S n (0,z) are equal 
to the modified approximants S n (w n , z) of the original continued fraction. The transformed 
continued fraction exists and can be calculated as 

d = b + w , Ci = Ai, d 1 = bi + w 1 , (4.13) 
Ci = Oi-i^-i, di = bi + Wi - Wi- 2 qi-i, i > 2, 
Xi = di - Wi-i(bi + Wi), qi = X i+1 /Xi i>l, 

if and only if A, ^ for i — 1,2, ... . 

We now return to the three-term recurrence relation focusing on their intimate relations 
to continued fractions. A three-term recurrence relation can be written as 

X n+l = b n X n + a n X n _! n = 1, 2, 3, . . . , (4.14) 

where complex numbers and a n ^ 0. The solutions of a three-term recurrence 

relation span a two-dimensional linear space. The {x n } nontrivial (i.e. ^ {0}) solution is 
said to be minimal if there exists another solution {y n } such that 

lim x n /y n = 0. (4-15) 

n^oo 

Solution {y n } here is called dominant. The minimal solution is unique, apart from a multi- 
plicative constant. 

The existence of the minimal solution is strongly related to the convergence of a continued 
fraction constructed from the coefficients of the recurrence relation. This connection is 
revealed by Pincherle's theorem [111,50]. According to this the following statements hold: 



A: ( [1.14j ) has a minimal solution if and only if the continued fraction 

K?- 1 (£) = £ £ £ (4.16) 
^KbnJ 6i + 6 2 H h 6„ H 

converges, 

B: if {X n } is a minimal solution then for AT = 0, 1, 2, ... , 

x n+i T ,oo ( a n +N\ ai+N a2+N a n+N 

= -A^i 7 = T . 7 . . 7 . • (4.17) 

XN \ n +N ) Ol+JV + »2+Af H ' + n+ N + " ' 

The second statement asserts that the ratio of two successive element of the minimal solution 
is provided by a continued fraction. 
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V. CONTINUED FRACTION FOR G 0N+ i/G, 



ON 



First we show that in certain domain of the complex plane the physical relevant solution 
of the recurrence relation ( |2.6| ) for the Green's matrix is the minimal solution. In case of 



short-range potentials the Green's function can be constructed as [22 



G{rS,k) = <p l (k,r < )fi + \k,r > )/F(k), (5.1) 

where (fi(k, r) is the regular solution, f\ + ' (k, r) is the Jost solution, ^(k) is the Jost function 
and k is the wave number. The Jost solution is defined by the relation 

lim e* ikr f\ ±] (k,r) = I. (5.2) 

Let us define a "new" Green's function as 

G(r, r', k) = <f t (k, r< )fi{k, r> )/F{k), (5.3) 

where fi is a linear combination of f\ + ^ and /[ '. If $t.E < /j is exponentially 
decreasing and f\ ^ is exponentially increasing. Thus, for any G we have 

lim G(r ' r/ ' fc) = 0, if m < 0. (5.4) 

r'-»oo G(r, r', k) 

We note, that both G and G satisfy the defining equation Eq. (|2.3j), but only G of Eq. ( |5.1p 
is the physical Green's function. The above considerations, with a slight modification in Eq. 
Q5.2|), are also valid for the Coulomb case. 

An interesting result of the study of Ref. |L3| is that the Green's matrix from Jacobi- 
matrix Hamiltonian, in correspondence with (|2.8|) , has an analogous structure to Eq. (|5.1| ) 

G w (k) = ((pi) i< (k)(fi + ^) i> (k)/ ^(k), (5.5) 

where J(pi)i(k) = ((pi(k)\i) and {fi +) )i{k) = (fj +) (k)\i) . Similarly, we define (fi)i(k) = 
(fi(k)\i) and 

G lll (k) = (< fl ) l< (k)(f l ) l> (k)/F(k). (5.6) 

On the dtE < region of the complex plane as r — > oo fi(k, r) exponentially dominates over 
fi + \k,r), thus for their L 2 representation the following relation holds 

lim ^tttt- = 0, if < 0. (5.7) 
{fi)i{k) 

This implies a similar relation for the Green's matrices 

l im = o, if m < 0. (5.8) 
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So, in the 3t.E < region of complex i?-plane the physical relevant Green's matrix Gu> 
appears as the minimal solution of recurrence relation ( |2.6|) . Thus, according to Pincherle's 
theorem ( 4.17 ), the ratio needed in Eq. ( |2.10| ) for the Green's matrix can be calculated by 
the continued fraction 



GW+i( e ) 
GW(e) 



-Kl 



--N 



(5.9) 



where dj = — Ju-i/ Ju+i, h = — Ju/Ju+i and J^- is the Jacobi-matrix. 
In the case of D-dimensional Coulomb Green's matrix we have 



i(i + 2l' + 1) 
\ (i + l)(i + 2Z' + 2) 



1,2, 



(5.10) 



m = n'-w+i i = 0>1|2> ... ( , n) 

(e + 6|) J(t + l)(t + 2l' + 2) 



with e = 2mE/h 2 , Z' = 2mZ/h 2 and t = I + (D - 3)/2. 

On the region of scattering states the recurrence relation does not have minimal solution 
and the continued fraction ( |5.9|) diverges. On the other hand, GW+i(e) / Gqn (e) is an analytic 
function of e, and there is a region of the complex plane where we have a representation for 
this function, thus values on other regions can be obtained by the analytic continuation of 
the bound-state formula, i.e. by the analytic continuation of the continued fraction (p.9|). 
Since now we have a limit 1-periodic continued fraction this, according to Eq. (|4.12| ), can 
be achieved with fixed points 

w ± (e) = -6(e)/2 ± y/ (b(e)) 2 / 4 + a (5.12) 

where 



a = lim a>i = — 1 (5.13) 

i— >oo 

6(e) = lim 6 i = 2(e-6|)/(e + 6|). 



Considering the formula for Green's operators f23 



(i\G(E + i0)\l) - (i\G(E - i0)\i) = -27ri(T|V>(£))(V>(£)K>, (5.14) 

where ip(E) is the scattering wave function, and the analytic properties of Green's operators 
we can readily derive that the imaginary part of (i\G(E + i0)|i) should be negative. This 
condition can only be fulfilled with the choice of w + . This choice gives an analytic contin- 
uation to the physical sheet, while u>_, which also converges, gives an analytic continuation 
to the unphysical sheet. 

^From the above considerations it follows that utilizing the Jacobi-matrix only the 
Green's matrix can be obtained for arbitrary complex energies by simply evaluating a con- 
tinued fraction. 
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VI. NUMERICAL ILLUSTRATIONS 



Below we demonstrate the convergence and the numerical accuracy of method B. We 
calculate the matrix element G^(e) of D-dimensional Coulomb Green's operator for the 
/ = and D = 3 case at bound- and scattering-state energies. We examine the convergence 
of continued fraction with different choice of {w n } in Eq. ( |4.6|) and the effect of Bauer-Muir 
transformation. For comparison we also give the exact value for G^ ( e ) (ELll)- 

In case of 3?e < we take the w n = 0, w n = w + and w n = W- choices. In Table I we can 
observe excellent convergence to the exact value in all cases. The choice of w n influences 
only the speed of convergence. 

In the region of 9fte > 0, in complete accordance with Pincherle's theorem, the contin- 
ued fraction ( |5.9| ) diverges, only its analytic continuation with w n = w + and w n = is 
convergent. However, as the first column in Table II shows, the convergence is rather poor. 
This can considerably be improved by the repeated application of Bauer-Muir transforms. 
In fact, an accuracy similar to the bound-state case can easily be reached here with e.g. an 
eightfold Bauer-Muir transform. 

In order that we have a more stringent test we have performed the contour integral 

1(C) = ^-i deG o(e). (6.1) 

If the domain surrounded by C does not contain any pole, then /(C) = 0. If this domain 
contains a single bound-state pole, then 1(C) = (0\ip) (ip\0) must hold, while if C circumvents 
the whole spectrum then 1(C) = (0 1 0) is expected. With appropriate selection of Gauss 
integration points we could reach 12 digits accuracy in all cases. This demonstrates that 
the calculation of the Green's matrices from J-matrices via continued fractions is accurate 
on the whole complex plane. 

VII. SUMMARY AND CONCLUSIONS 

In this paper we have shown that if in some basis representation the Hamiltonian takes a 
Jacobi-matrix form the corresponding Green's matrix can be calculated on the whole com- 
plex energy plane by a continued fraction, whose coefficients are related to the elements of 
the Jacobi-matrix. To justify this statement we presented the example of the D-dimensional 
Coulomb problem, in particular, we calculated the Coulomb-Sturmian-space representation 
of the D-dimensional Coulomb Green's operator. Numerical examples proved the accuracy 
and the efficiency of the method. 

The applicability of the techniques presented here can be extended beyond the examples 
discussed above. We may have a physical situation in which only the asymptotic part of 
the Hamiltonian has Jacobi-matrix structure. In this case, like in the L 2 approach we can 
approximate V by its projection onto a finite subspace spanned by the first N basis 
states and can generate the analytic Green's matrix as a solution of a Lippmann-Schwinger 
matrix equation (see e.g. in Ref. MM)- Also, our analytic two-body Green's matrices may 
be used to derive Green's matrices of composite systems via convolution integrals based on 
two-body problems. 
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It should be emphasized that we used the Jacobi-matrix form of a particular Hamil- 
tonian only and the method is applicable to any Jacobi-matrix Hamiltonian if the matrix 
elements are know analytically. This later requirement may be relaxed and thus we can de- 
termine approximate Green's matrices from J-matrices generated by the numerical Lanczos 
procedure. 
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VIII. APPENDIX 

Here we present the formulas analogous to those in Section [TTI] for the D-dimensional 
Harmonic oscillator. The radial Schrodinger equation in this case is 



H H °^{r) 
The energy eigenvalues are 



2m \ dr 2 r 2 V 2 A 2 J 2 



ip(r) = Eip(r). (8.1) 



E nl = hu(2n + l + -\ (8.2) 



2 

and the corresponding wave functions can be written as 

1 /2 

/i ,v , / x 1 / 2T(n+l) \ , v 9 , , 9 . i , d^x , n+R -i) , 9n 

(r u,nl) = i[> nl (uj,r) = v-*[ \ J D . ) exp(--r 2 )(vr 2 ) LV 2 \vr% (8-3) 

\T(n + l + ^)J 2 

where v = mu/Ti. The harmonic oscillator functions are orthonormal and form a complete 
set in the usual sense. 

The harmonic oscillator Hamiltonian with parameter uo on the basis of harmonic oscillator 
function with parameter u' takes a Jacobi-matrix form 

<u/, nl\H HO (uj, E) |u/, n'l) = 5 nn , ^ fa + 1 + ^ 

. X-^' 2 ( ,( , , ' D A\ 1/2 

-5 nn/+1 h^-^ fa + l)(n' + l+ ^)) V2 (8.4) 

The calculation of the Green's matrix via method B goes analogously to the Coulomb 
case. The matrix element (0l\G HO (E)\0l), which can be used in method A, is given as 

G ,o (£)|0iH _ 8 ^_J_ 

(iD E I D E (uj-uj'\ 2 \ , . 

x 2 F a + i_ I- _ + _ + !_ • . 8.5 

2 1 \ 2 4 2hu '2 4 2ftw' \u + uj' 
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TABLES 



n 


c4 0) 


e = (-100, 0) 




1 


(-5.44922314793965,0) 


(-5.59142801316938,0) 


(-0.92906408986331,0) 


2 


(-5.54075476366523,0) 


(-5.56131039101044,0) 


(-4-70080363351349,0) 


3 


(-5.55501552420656,0) 


(-5.55812941271530,0) 


(-5.39340957492282,0) 


4 


(-5.55726787304507,0) 


(-5.55775017508704,0) 


(-5.52662100417805,0) 


5 


(-5.55762610832912,0) 


(-5.55770176067796,0) 


(-5.55192403535264,0) 


6 


(-5.55768333962797,0) 


(-5.55769530213874,0) 


(-5.55663951922343,0) 


7 


(-5.55769251168083,0) 


(-5.55769441374319,0) 


(-5.55750389878413,0) 


8 


(-5.55769398510141,0) 


(-5.55769428874846,0) 


(-5.55766025755765,0) 


9 


(-5.55769422223276,0) 


(-5.55769427085427,0) 


(-5.55768824220786,0) 


10 


(-5.55769426045319,0) 


(-5.55769426825710,0) 


(-5.55769320762502,0) 


11 


(-5.55769426662100,0) 


(-5.55769426787592,0) 


(-5.55769408236034,0) 


12 


(-5.55769426761735,0) 


(-5.55769426781946,0) 


(-5.55769423553196,0) 


13 


(-5.55769426777843,0) 


(-5.55769426781103,0) 


(-5.55769426221577,0) 


14 


(-5.55769426780450,0) 


(-5.55769426780976,0) 


(-5.55769426684377,0) 


15 


(-5.55769426780872,0) 


(-5.55769426780957,0) 


(-5.55769426764335,0) 


16 


(-5.55769426780940,0) 


(-5.55769426780954,0) 


(-5.55769426778103,0) 




(-5.55769426780951,0) 


(-5.55769426780954,0) 


(-5.55769426780465,0) 


18 


(-5.55769426780954,0) 




(-5.55769426780870,0) 


19 


(-5.55769426780954,0) 




(-5.55769426780939,0) 


20 






(-5.55769426780950,0) 


21 






(-5.55769426780954,0) 


22 






(-5.55769426780954,0) 




G A = (-5.55769426780954, 0) 



TABLE I. Convergence of the continued fraction for first element of the Green's matrix at 
!fte < with method B. The first, second and third column contain approximants of the continued 
fraction with w n = 0, w n = w+ and w n = W-, respectively. For comparison we also give the exact 
result. All the G values are scaled with 10 2 . 
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(w+) 



G B 
(1.076. 
(1.074 
(1.198 
(1.110 
(1.160 
(1.141 
(1.139 
(1.157 
(1.131 
(1.158 
(1.135 
(1.149 
(1.145 
(1.140 
(1.152 
(1.137 
(1.152 
(1.140 
(1.147 
(1.146 



(0) 
■0.678) 
•0.279) 
•0.325) 
■0.346) 
■0.307) 
■0.354) 
■0.312) 
■0.341) 
■0.325) 
■0.327) 
■0.337) 
■0.320) 
■0.340) 
■0.322) 
■0.334) 
■0.329) 
■0.327) 
■0.335) 
■0.323) 
■0.336) 



(1.8072 
(1.1225 
(1.1425 
(1.1497 
(1.1415 
(1.1478 
(1.1437 
(1.1458 
(1.1451 
(1.1448 
(1.1457 
(1.1446 
(1.1456 
(1.1449 
(1.1453 
(1.1452 
(1.1450 
(1.1454 
(1.1450 
(1.1453 



0>+), 



e = (1000, 1) 



1) 



(w+) 



,-0.3293) 
,-0.3783) 
,-0.3162) 
,-0.3353) 
,-0.3287) 
,-0.3298) 
, 0.3311) 
,-0.3290) 
,-0.3312) 
,-0.3294) 
,-0.3305) 
,-0.3300) 
,-0.3300) 
,-0.3304) 
,-0.3298) 
,-0.3304) 
,-0.3296) 
,-0.3302) 
,0.3301) 
,-0.3300) 



G_B_ 

(-0.4129544,- 

(4.29352799, 
(1.13445656, 
(1.14598003, 
(1.14512823, 
(1.14523860, 
(1.14522511, 
(1.14522377, 
(1.14522642, 
(1.14522458, 
(1.14522559, 
(1.14522512, 
(1.14522528, 
(1.14522527, 
(1.14522523, 
(1.14522528, 
(1.14522524, 
(1.14522527, 
(1.14522525, 
(1.14522526, 



J5) 

■0.14238595) 
-1.63424931) 
-0.32244006) 
-0.33019962) 
-0.33023791) 
-0.33015825) 
-0.33018993) 
-0.33017859) 
-0.33018226) 
-0.33018138) 
-0.33018133) 
-0.33018161) 
-0.33018135) 
-0.33018153) 
-0.33018143) 
-0.33018147) 
-0.33018146) 
-0.33018145) 
-0.33018147) 
-0.33018145) 



G 



(w+) 
B 



(-0.2321154,-0 
(-1.4408861,-0 
(1.20667672,0. 
(1.14702383,-0 
(1.14511731,-0 
(1.14523597,-0 
(1.14522395,-0 
(1.14522539,-0 
(1.14522527,-0 
(1.14522524,-0 
(1.14522527,-0 
(1.14522525,-0 
(1.14522526,-0 
(1.14522525,-0 
(1.14522526,-0 
(1.14522526,-0 
(1.14522526,-0 
(1.14522526,-0 
(1.14522526,-0 
(1.14522526,-0 



(8) 

.073120618) 
.350899497) 
237375310) 
.332562329) 
.330140243) 
.330179581) 
330182552) 
330181124) 
330181563) 
330181440) 
330181470) 
.330181465) 
.330181464) 
.330181466) 
.330181465) 
.330181466) 
.330181465) 
.330181465) 
.330181465) 
330181465) 



G A = (1.14522526, -0.330181465) 



TABLE II. Convergence of the continued fraction for first element of the Green's matrix at 
Jfte > with method B. The first, second, third and fourth column contain approximants of the 
continued fraction with w n = w + and without, with one-fold, with five-fold and with eight-fold 
Bauer-Muir transform, respectively. For comparison we also give the exact result. All the Goo 
values are scaled with 10 2 . 



16 



